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Abstract 

Building quantum devices using fixed operators is a must to simplify the hardware con- 
struction. Quantum search engine is not an exception. In this paper, a fixed phase quantum 
search algorithm that searches for M matches in an unstructured search space of size N 
will be presented. Selecting phase shifts of 1.916847T in the standard amplitude amplification 
will make the technique perform better so as to get probability of success at least 99.58% in 




better than any know fixed operator quantum search algorithms. The algorithm 



will be able to handle either a single match or multiple matches in the search space. The 
algorithm will find a match in O \ ^J~N]~M\ whether the number of matches is known or not 
in advance. 

1 Introduction 

In 1996, Lov Grover [10] presented an algorithm that quantum mechanically searches an unstruc- 
tured list assuming that a unique match exists in the list with quadratic speed-up over classical 
algorithms. To be able to define the target problem of this paper, we have to organize the ef- 
forts done by others in that field. The unstructured search problem targeted by Grover's original 
algorithm is deviated in the literature to the following four major problems: 

• Unstructured list with a unique match. 

• Unstructured list with one or more matches, where the number of matches is known 

• Unstructured list with one or more matches, where the number of matches is unknown. 

• Unstructured list with strictly multiple matches. 

The efforts done in all the above cases, similar to Grover's original work, used quantum paral- 
lelism by preparing superposition that represents all the items in the list. The superposition could 
be uniform or arbitrary. The techniques used in most of the cases to amplify the amplitude (s) of 
the required state(s) have been generalized to an amplitude amplification technique that iterates 
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the operation UR S (0) U^R t (<p), on U \ s) where U is unitary operator, R s (4>) = I — (1 — e % ^) \s) (s\, 
R t (<f) = I — (1 — e tip ) \t) (t\, \s) is the initial state of the system, \t) represents the target state(s) 
and / is the identity operator. 

Grover's original algorithm replaces U be W, where W is the Walsh-Hadamard transform, pre- 
pares the superposition W |0) (uniform superposition) and iterates WR S (ir) WR t (ir) for O ( Vn) , 



where N is the size of the list, which was shown be optimal to get the highest probability with the 
minimum number of iterations [23], such that there is only one match in the search space. 

In [TT1 [T5l [HI [17j [l], Grover's algorithm is generalized by showing that U can be replaced by 
almost any arbitrary superposition and the phase shifts <j) and <p can be generalized to deal with the 
arbitrary superposition and/or to increase the probability of success even with a factor increase 
in the number of iterations to still run in 0(\/N). These give a larger class of algorithms for 
amplitude amplification using variable operators from which Grover's algorithm was shown to be 
a special case. 

In another direction, work has been done trying to generalize Grover's algorithm with a uniform 
superposition for known number of multiple matches in the search space [31 [HI Q E] , where it was 
shown that the required number of iterations is approximately 7r/4^/ N/M for small M/N, where 
M is the number of matches. The required number of iterations will increase for M > N/2, 
i.e. the problem will be harder where it might be excepted to be easier [19]. Another work has 
been done for known number of multiple matches with arbitrary superposition and phase shifts 
[T5J [21 HI dH [TH] where the same problem for multiple matches occurs. In [3 dBJ H], a hybrid 
algorithm was presented to deal with this problem by applying Grover's fixed operators algorithm 
for w/&WN~/M times then apply one more step using specific and cp according to the knowledge 
of the number of matches M to get the solution with probability close to certainty. Using this 
algorithm will increase the hardware cost since we have to build one more R s and R t for each 
particular M . For the sake of practicality, the operators should be fixed for any given M and are 
able to handle the problem with high probability whether or not M is known in advance. In [2U [22] , 
Younes et al presented an algorithm that exploits entanglement and partial diffusion operator to 
perform the search and can perform in case of either a single match or multiple matches where 
the number of matches is known or not [22] covering the whole possible range, i.e. 1 < M < N. 
Grover described this algorithm as the best quantum search algorithm [12]. It can be shown that 
we can get the same probability of success of [21] using amplitude amplification with phase shifts 
<fi = if = it/2, although the amplitude amplification mechanism will be different. The mechanism 
used to manipulate the amplitudes could be useful in many applications, for example, superposition 
preparation and error-correction. 

For unknown number of matches, an algorithm for estimating the number of matches (quantum 
counting algorithm) was presented [SKIS]- In [3], another algorithm was presented to find a match 
even if the number of matches is unknown which will be able to work if M lies within the range 
1 < M < 3N/A [22]- 

For strictly multiple matches, Younes et al [20] presented an algorithm which works very ef- 
ficiently only in case of multiple matches within the search space that splits the solution states 
over more states, inverts the sign of half of them (phase shift of -1) and keeps the other half 
unchanged every iteration. This will keep the mean of the amplitudes to a minimum for multiple 
matches. The same result was rediscovered by Grover using amplitude amplification with phase 
shifts (f) = f = 7r/3 [13], in both algorithms the behavior will be similar to the classical algorithms 
in the worst case. 



This algorithm is able to handle the range 1 < M < N for both known and unknown number of 




In this paper, we will propose a fixed phase quantum search algorithm that runs in O 
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matches more reliably than known fixed operator quantum search algorithms that target this case. 

The plan of the paper is as follows: Section 2 introduces the general definition of the target 
unstructured search problem. Section 3 presents the algorithm for both known and unknown 
number of matches. The paper will end up with a general conclusion in Section 4. 

2 Unstructured Search Problem 

Consider an unstructured list L of N items. For simplicity and without loss of generality we will 
assume that iV = 2" for some positive integer n. Suppose the items in the list are labeled with the 
integers {0, 1, ...,N — 1}, and consider a function (oracle) / which maps an item i G L to either 
or 1 according to some properties this item should satisfy, i.e. / : L — > {0, 1}. The problem 
is to find any i 6 I such that f(i) = 1 assuming that such % exists in the list. In conventional 
computers, solving this problem needs O (N/M) calls to the oracle (query), where M is the number 
of items that satisfy the oracle. 

3 Fixed Phase Algorithm 
3.1 Known Number of Matches 

Assume that the system is initially in state \s) = |0). Assume that ^/ denotes a sum over % 
which are desired matches, and denotes a sum over % which are undesired items in the list. 
So, Applying U \ s) we get, 

N-l , N-l 

|</>< 0, > = c| s > = ^£> : > + ^£"l°>. ID 

V iV i=0 v i=0 

where U = W and the superscript in represents the iteration number. 

Let M be the number of matches, sin(6 l ) = \J M/N and < 6 < tt/2, then the system can be 
re-written as follows, 

|^ (0) > =sm(0)|V>i>+cos(0)|^ o >, (2) 

where \ipi) = \t) represents the matches subspace and |^o) represents the non-matches subspace. 

Let D = UR S (0) WR t (<p), R s (<p) = I - (1 - e**) \s) (s\, R t (<p) = I - (1 - e*) \t) (t\ and set 
<p = (p as the best choice [H]. Applying D on we get, 

|^ 1 )> = J D|^ (0) > = a 1 |Vi}+6i|^o), (3) 

such that, 

ai=sin(0)(2cos(<J)e^ + l), (4) 

6i = e < *cos(0)(2cos( < J) + l), (5) 

where cos (5) = 2 sin 2 (6) sin 2 (|) — 1. 

Let q represents the required number of iterations to get a match with the highest possible 
probability. After q applications of D on |?/^ ^ we get, 

|^> = ^|^ 0) > = a g ^i> + & 9 |^o>, (6) 



3 



such that, 



a q = sin 



(9)(e^U q (y) + e^- 1 ^U q ^(y)), (7) 



b q = cos(0)e^-^ (U q (y) + U q -! (y)) , (8) 
where y = cos(5) and U q is the Chebyshev polynomial of the second kind defined as follows, 

sin ((g +l)g) 

" ,U/J sin (5) ' W 

Let P 9 represents the probability of success to get a match after q iterations and P 9 S is the 

i 1 2 i 1 2 

probability not to get a match after applying measurement, so P 9 = \a q \ and P 9 S = \b q \ such 
that P 9 + P£ s = 1. To calculate the required number of iterations q to get a match with certainty, 
one the following two approaches might be followed: 

• Analytically. The usual approach used in the literature when the number of matches M is 
known in advance is to equate P 9 to 1 or P 9 S to and then find an algebraic formula that 
represents the required number of iterations, as well as, the phase shifts and (p in terms on 
M. Using this approach is not possible for the case that the phase shifts should be fixed for 
an arbitrary M such that 1 < M < N as shown in the following theorem. 

Theorem 3.1 (No Certainty Principle) Let D be an amplitude amplification operator 
such that D = UR S (0) U^R t ((f), where U is unitary operator, R s (<p) = I — (1 — e 1 ^) \s) (s\, 
R t (f) = I — (1 — e l(p ) \t) (t\, \s) is the initial state of the system, \t) represents the target 
state(s) and I is the identity operator. Let D performs on a system initially set to U\s). 
If the phase shifts <fr and (p should be fixed, then iterating D an arbitrary number of times 
will not find a match with certainty for an arbitrary known number of matches M such that 
1 < M < N. 

Proof To prove this theorem, we will use the usual approach, i.e. start with P 9 = 1 or 
P% s = and calculate the required number of iterations q. 

Since Pf = \a q \ 2 and from Eqn^fy we can re-write P s 9 as follows setting <f> — if as the best 
choice [L$ , 

™ sin 2 (9) , 



sin 2 (S) 



cos (5) cos ((2q + 1) 5) + 2 cos (0) sin ((q + 1) 5) sin (qS)) . (10) 



Setting P 9 = 1 and using simple trigonometric identities we get, q = i.e. the required 
number of iterations is independent of M, <p and tp, and represents an impossible value for a 
required number of iterations. 

Direct Search. The alternative approach used in this paper is to empirically assume an 
algebraic form for the required number of iterations that satisfy the quadratic speed-up of 
the known quantum search algorithms and use a computer program to search for the best 
phase shift (f> that satisfy the condition, 

max (min (P s 9 (0))) such that < < 2ir and, 1 < M < N. (11) 
i.e. find the value of <fi that maximize the minimum value of P 9 over the range 1 < M < N . 
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Figure 1: The probability of success the proposed algorithm after the required number of iterations. 



written using C language to find the best with precision 10~ 15 that satisfy the conditions shown 
in Eqn. [TT] The program shows that using = 6.021930660106538 ~ 1.916847T, the minimum 
probability of success will be at least 99.58% compared with 87.88 % for Younes et al [22] and 
50% for the original Grover's algorithm [3] as shown in Fig. [2j To prove these results, using 
= 1.916847T, the lower bound for the probability of success is as follows as shown in Fig. [TJ 



P s = S$ (1 - cos W cos ^ 2q + !) 5 ) + 2 cos W sin ((<? + !) 6 ) sin (?<*)) 
= (1 - cos (5) cos ((2q + 1) 5) + cos (0) cos (5) - cos ((f)) cos ((2q + 1) 5)) (12) 

> (i + cos 2 ( 5 ) + 2 cos (0) cos (5)) > 0.9958. 

where, cos (5) = 2sin 2 (fl) sin 2 (f ) - 1, < 9 < vr/2, and cos ({2q + 1)S)< -cos(5). 
3.2 Unknown Number of Matches 



In case we do not know the number of matches M in advance, we can apply the algorithm shown 
in [3] for 1 < M < N by replacing Grover's step with the proposed algorithm. The algorithm can 
be summarized as follows, 

1- Initialize m = 1 and A = 8/7. (where A can take any value between 1 and 4/3) 

2- Pick an integer j between and m — 1 in a uniform random manner. 

3- Run j iterations of the proposed algorithm on the state |?/;( ^: 



Assume that q 



sin(6») 




Using this form for g, a computer program has been 



^0')) = D j |^ (0) ) . 



(13) 



1.2 - 




0.4 - 



0.2 











0.2 



0.4 



0.6 



0.8 



M/N 



Figure 2: The probability of success of Grover's algorithm, Younes et al algorithm [21] and the 
proposed algorithm after the required number of iterations. 



where m represents the range of random numbers (step 2), j represents the random number of 
iterations (step3), and A is a factor used to increase the range of random numbers after each trial 
(step 6). 

For the sake of simplicity and to be able to compare the performance of this algorithm with 
that shown in [3], we will try to follow the same style of analysis used in [3]. Before we construct 
the analysis, we need the following lemmas. 

Lemma 3.2 For any positive integer m and real numbers 6, 5 such that cos (5) = csin 2 (6 l ) — 1, 
< 6 < 7r/2 where c = 2sin 2 (|) is a constant, 
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Measure the register W)^>} and assume i is the output. 
If f{i) = 1, then we found a solution and exit. 





q=0 



Proof By mathematical induction. 



Lemma 3.3 For any positive integer m and real numbers 6 
< 6 < it/2 where c = 2sin 2 (|) is a constant, 



5 such that cos (5) = csin 2 (6') — 1 




(3 



Proof By mathematical induction. 



Lemma 3.4 Assume M is the unknown number of matches such that 1 < M < N. Let 9, 5 be real 
numbers such that cos (5) = 2 sm 2 (6) sin 2 (f ) - I, sm 2 (6) = M/N, <p = 1.91684vr and < 9 < ix/2. 
Let m be any positive integer. Let q be any integer picked in a uniform random manner between 
and m — 1. Measuring the register after applying q iterations of the proposed algorithm starting 
from the initial state, the probability P m of finding a solution is as follows, 



1 / ... . .. (cos (5) + cos ((f))) sin (2m8)\ 

— 1 + cos (5) cos ((f)) - i ^— . VV 7; '- 

cos(oj) \ 2msm{d) J 



c (1 — cos (5)) 

where c = 2sin 2 (|) ; then P m > 1/4 form > 1/ sin (5) and small M/N. 

Proof The average probability of success when applying q iterations of the proposed algorithm 
when < q < m is picked in a uniform random manner is as follows, 

m—1 

P = — V P q 

q=0 

n , . m—1 
sin 2 (6>) 



m sir?(5) 

q=0 



]T (sin 2 ((q + 1)5) + sin 2 (qS) + 2 cos ((f)) sin ((q + 1) 5) sin (qS)) 



sin 2 (9) f _ cos((5) sm{2m5) / ,\ / « _ cos (0) sin(2m<S) \ 

msin 2 (5) y a 2sin(«5) ~r COS \(p) COS yu ) 2sin(5) J 



± ' 1 + COS (6) COS ((f)) - (c°sffl+c°sW)sin( 2m ,S) 



msin (5) 
1 

c{\— cos(<5)) V 1 VjVJ,J V" i Vj>JlJ \f l 2msin(i5) 

If m >!/ sin (S) and M -C N then cos (5) ~ —1, so, 
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If . . (cos (6)- 1) sin (2mo1 \ 1 / . . (l-cos(0))\ 
P m >_h -cosfo)-* ^ >-) >-(l-cos(0)-^ =0.25 

where — 1 < sin (2mo~) < 1 /or < # < 7r/2. 

We calculate the total expected number of iterations as done in Theorem 3 in [3] . Assume that 
m q > 1 1 sin (5), and v q = [log A m q ] . Notice that, m q = O \^/N/Mj for 1 < M < N, then: 

1- The total expected number of iterations to reach the critical stage, i.e. when m > m q \ 

1 v " i 

- V A'"" 1 < — -m q = 3.5m„. (14) 

2 ^ ~ 2 (A - 1) 9 q V ' 

2- The total expected number of iterations after reaching the critical stage: 



1 00 /Q\ n 

;Ed a" ,+u — ______ my — 3.5 m? . (1B) 



2^14/ 2(1-0.75A) 

The total expected number of iterations whether we reach to the critical stage or not is 7m q 
which is in 0(y/N/M) for 1 < M < N. 

When this algorithm employed Grover's algorithm, and based on the condition mc > 1/ sin (29 g) 

yy/N/Mj for M < 3iV/4,the total expected number of iterations is approximately 8mc for 

1 < M < 3N/4. Employing the proposed algorithm instead, and based on the condition 
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Figure 3: The actual behavior of the functions representing the total expected number of iterations 
for Grover's algorithm, Younes et al algorithm [22] and the proposed algorithm taking A = 8/7, 
where the number of iterations is the flooring of the values (step function). 



m q > 1/ sin (5) = O y^/N/MJ ,the total expected number of iterations is approximately 7m q 

for 1 < M < N, i.e. the algorithm will be able to handle the whole range, since m q will be able to 
act as a lower bound for q over 1 < M < N . Fig. [3] compares between the total expected number 
of iterations for Grover's algorithm, Younes et al algorithm [22] and the Fixed Phase algorithm 
taking A = 8/7. 



4 Conclusion 

To be able to build a practical search engine, the engine should be constructed from fixed operators 
that can handle the whole possible range of the search problem, i.e. whether a single match or 
multiple matches exist in the search space. It should also be able to handle the case where the 
number of matches is unknown. The engine should perform with the highest possible probability 
after performing the required number of iterations. 

In this paper, a fixed phase quantum search algorithm is presented. It was shown that selecting 
the phase shifts to 1.916847T could enhance the searching process so as to get a solution with 
probability at least 99.58%. The algorithm still achieves the quadratic speed up of Grover's 
original algorithm. It was shown that Younes et al algorithm [22] might perform better in case 
the number of matches is unknown, although the presented algorithm might scale similar with an 

acceptable delay, i.e. both run in O (yN/M\ In that sense, the Fixed Phase algorithm can act 

efficiently in all the possible classes of the unstructured search problem. 
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